function[quiet_jicamarca_response] = jic_quiet(fday)
% The function returns the average, quiet day response of the
% vertical drift at Jicamarca Radar site
% fday should be in decimal day, in UT

if ispc == 1,
    
    load E:\projects\ace_tensor\jicamarca_quiet_average sp;
    
end;

if isunix == 1,
    
    load /home/mnair/projects/ace_tensor/jicamarca_quiet_average.mat sp;
    %   eval(['load ' data_dir 'ace_tensor/jicamarca_quiet_average sp']);
       
end
    

fday_decimals_ = fday - floor(fday);

quiet_jicamarca_response = ppval(fday_decimals_,sp);
